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The propagation of audiofrecjuency and sub-audiofrequency waves between the earth 
and an ionosphere whose conductivity varies continuously with altitude is considered in 
detail. The fields are represented in terms of two scalar potentials satisfying appropriate 
w^ave equations in spherical coordinates. 

It is shown, on the basis of existing data on the ionosphere, that waves in this frequency 
range can l)e considered to be confinecl to a thin, but not sharply Ixninded, spherical shell 
about the earth. Greatly simplified radial wave equations in dimensionless form are 
deriv^ed incorporating this approximation. Solutions of these equations are given for two 
regions, viz, for the low-altitude region where (T/coeo<<l and, in the case of certain 
restricted types of conductivity profile, for the high-altitude region where o-/coe()>>l. 
An iterative method is presented, based on an integral (Hjuation, which makes possible 
a computation of the radial wave function in the transition region and a joining of interior 
and exterior solutions for the ])ropagating TM mode. The result is a direct nuithematical 
relationship between the conductivity ])rofile and the complex propagation constant as 
a function of frequency. 

It is demonstrated that at frecjuencies abo\-e about 50 cycles the ])ropagalion constant 
can be obtained from measurements of the horizontal components of electric and magnetic 
fields in individual sf erics at airplane altitudes, while at lower frequencies the same infor- 
mation can be obtained through ground-based observation of cavity resonance effects in 
sferic noise. Existing data on these cavity resonances are used to calculate the complex 
propagation constant as a function of frecpiency from 6 to 34 cycles. The results, when 
extra])olated to higher frequencies, predict attenuation rates in excellent agreement with 
currently available data. 

The effect of the diurnal variation in conductivity on observable* quantities is l)riefly 
examined, and tentative conclusions as to its magnitude are drawn. 

1. Introduction 

The propagation characteristics of ELF and sub-audiofrequency radio waves in the 
earth-ionosphere cavity are critically dependent upon the detailed electrical properties of the 
ionosphere. Most previous theoretical treatments [Budden, 1957; Wait, 1957; Wait, 1058; 
Wait, 1960 a and b] of the problem have assmned an ionosphere sliarply bomuled at some 
definite height, with a conductivity wliich, from that height outward, either remained constant 
or varied stepwise or continuously with height. For a number of reasons, such models can 
never adequately represent the actual physical situation. The same observation has been 
made recently by Wait [1962] in a paper treating the case of a continuously varying conduc- 
tivity which becomes constant at high altitude. The work presented below, while confining 
itself substantially to the ELF region, is an attempt to approach a somewhat more general 
problem from a different point of view. In particular, this paper is directed toward the 
possibility of inferring the properties of the lower ionosphere from field measurements taken 
over a range of frequencies. 

For frequencies low compared to the electronic collision frequency — a condition easily 
satisfied in the VLF range and below — the ionosphere behaves like a medium with unit 
dielectric constant, imit relative permeability, and a conductivity which is a function of 
position. The conductivity affects the local behavior of the electromagnetic field only as 
it enters into the complex relative permittivity e, which can be represented as follows: 



• Tills work was sponsotvd in part by the Geophysical Research Directorate, Air Force Cambridge Research Laboratories, Bedford, Mass. 
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e=l-i^' (1) 
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where cor is a real number, related to the conductivity, which depends upon position but not 
frequency. At altitudes where (cor/co) is much smaller than unity, the medium behaves 
essential^ like a vacuum. Where this ratio is much larger than unity, it behaves like a good 
conductor. Although (1) makes no provision for the effect of the terrestrial magnetic field, 
it has been shown [Wait, 1960 b] that this can be quite simply taken into account at low 
frequencies through the use of an effective co,-, dependent on the vertical component of magnetic 
field. 

The transition region, in which (co;./co) is comparable with or at most a few orders of magni- 
tude greater than unity, can play as significant a role in determining the propagation of ELF 
waves as does the upper region. It is precisely this region which is omitted from sharply 
bounded models. Moreover, it can be seen from the form of (1) that the height at which the 
effect of finite conductivity is first felt is progressively lower for waves of lower frequencA^ 
On the other hand, because of a phenomenon analogous to the familiar skin effect, the total 
penetration of the waves into the ionosphere increases with decreasing frequency. When 
viewed in this light, the concept of a single ''ionospheric height,^' even if this height is per- 
mitted to vary with frequency, clearly loses its physical meaning and probably its mathe- 
matical utihty. Certainly some effort to take into account a continuous variation of oor with 
height is warranted. 

Any theory which attempts to relate the conductivity profile of the ionosphere to observed 
fields must take into account both the real and the imaginary parts of the effective propagation 
constant as functions of frequency. Unfortunately experimental data [Chapman and Macario, 
1956; Jean, 1961] to date, with few exceptions [Balser and Wagner, 1960a and b], are capable 
of giving information only about the attenuation characteristics of the wave and yield nothing 
about its phase velocity. Future experiments on ELF and sub-audiofrequency radiation, if 
they are to lead to increased knowledge of ionospheric properties, must provide for measure- 
ment of both attenuation and phase lag. 

2. Scalar Solutions for the Spherically Symmetrical Ionosphere 

In the oversimplified, but important, case in which the conductivity depends only on the 

radius of spherical coordinates, solutions of Maxwell's equations can be found in the form of 

-> -> 

TM waves (r-component of H absent) and TE waves (r-component of E absent). Each type 

of solution can be derived from a single scalar potential. For the TM wave, the fields are 

given by 

H=rXVyP (2) 

E=y^ (3) 

where i/' is a scalar solution of 

^'^•(^^V')+/:V=0. (4) 

The symbol k is here used to represent the free-space propagation constant of a wave of the 
given frequency. For the TE wave, one obtains 

E=rXVx (5) 

^=-^ (6) 

where x is a scalar solution of 
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7^-(;|vx)+Fx=o. (7) 

Both (4) and (7) are separable in spherical coordinates, giving solutions of tlie lorm 

\^=z:i:<*,«we?wi?.(r) (8) 

v m 
At m 

The functions ^mW are linear combinations of exp {±jm(l>) where m is zero or an integer, 
and the G^orM(^) ^^^ appropriate Legendre functions. The radial parts of the solutions satisfy 
the following equations: 






-[f.-'^]s.=0. (11) 



The permissible values, or eigenvalues, v and n for a given frequency are determined b\' 
two boundary conditions. These are, in general, that tlie fields remain finite at the center of 
the earth and vanisli at infinity. In the case of a perfectly conducting earth, the first condi- 
tion is replaced by 



/dRA 

\drjr 



S,{ro) = (13) 

where /q is the radius of the earth. It should be noted that the eigenvalues for a pailicnhir 
frequency are uniquely determined by the radial equation alone, and are in gcMuvral ('()nij)]ex. 
The Legendre functions, which can be written PTor^^O: ^*os 6), will tlierefore in general possess 
at least one singularity (at d=T if tlie + sign is used, or ^=0 if the — sign is used), correspond- 
ing to the fact that a source is needed to drive the system at an arbitrary frequency. There 
will, however, be certain (complex) frequencies at which one of the eigenvalues becomes equal 
to an integer n. Since the PJJ^(cos 6) liave no singularities, these are the natural resonant fre- 
quencies of the system, at which oscillation is possible in the absence of a source. 

3. Altitude Dependence of the Relative Permittivity 

The parameter co^ depends on the electron density N and the collision frequency Vc as 
follows : 



fom Kvc / 



(14) 



where e and m are the charge and mass of the electron. Existing data [Pierce, 1960; Ratcliffe. 
1959] on A^ and Vc have been used to construct the curve of figure 1, which shows the probable 
approximate variation of a;,, with height above the surface of the earth during daylight hours. 
It is seen that the function, although continuous, is very rapidly varying and reaches extremel}^ 
higli values in the first 100 km. The increase continues thereafter, rather more slowly, until 
a broad peak at an altitude of 300 or 400 km is reached, where cor is several orders of magnitude 
larger than at 100 km. A relatively slow decrease then sets in, which results in a significant 
reduction of co^ onl}^ after an altitude of several thousand kilometers is reached. 
If the conductivity, given by 

a=€o(JOr (15) 

were even to remain constant from 100 km upward for several hundred kilometers, the skin 
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Figure 1. Approximate daytime ionosphere 
profile. 



depth in that region for ELF radiation would be very small. At the very low frequency of 
10 cycles, the skin depth would be only of the order of 1 km. Therefore, it would seem safe to 
assume that the wave is attenuated essentially to zero within a vertical distance which is small 
compared with the radius of the earth. Although figure 1 may be in error by an order of 
magnitude or more at a given altitude, and although the nighttime profile is somewhat different, 
certainly one can still be confident that all significant phenomena are confined to a shell not 
more than 100 to 150 km in thickness. 

The argument just outlined, together with the fact that the attenuation is also large in the 
earth itself, leads to an important mathematical simplification in the radial wave equations, 
which may now be written 









(16) 
(17) 



The variable r in the final terms has been replaced by a constant effective radius f^ which can 
be chosen to be somewhat larger than Vq. The form of the relevant solutions will be practically 
unaffected by this change, even if the independent variable is allowed to go to infinity in 
satisfying the outer boundary condition. 

The simplified forms of the radial equations make it easy to attach a direct ph^^sical mean- 
ing to the eigenvalues. Equations (16) and (17) are identical to those which are obtained for 
the potentials in the case of a wave propagating between infinite plane-parallel w^alls having 
the same electrical characteristics as those of the earth and the ionosphere, providing only that 
the following replacements are made: 

=2 "^^TM (18) 
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'^=)tk (19) 

where k^yi and /:te iH'o the propagation constants of TAI and TE waves, respcH'tivc^ly, in a 
direction parallel to the walls. 

4. Dimensionless Form of the Radial Equations 

For simplicity of notation, it is convenient to rewrite (18) and (19) in dimensionless form, 
as follows: 

^+(«-/3)r,=0 (21) 

where 

x=.k{r-n) (22) 

— ^^^ + ^^ (23) 



Ua(x)=h\{r) (25) 

v,{x)=S,{r). (26) 

In order to see approximately the nature of the solutions and the order of magnitude of the 
eigenvalues, it is interestuig to consich^r the case of perfectly con(hicting earth ami ionosphere, 
separated by a vacuum of thickness h. The solutions would then he of the foi'ni 

Ua=cos Wl—ax] (27) 

?'^=sin yi-^xl (28) 



The bouTulary conditions at x=kh reciuire that 



Vl-a=g (p=0,l,2,3. . .) (29) 



•V' 



_qir 



/l-^=g (2=1,2,3...). (30) 



Since, at frequencies of the order of those under discussion, kh is much smaller than 
unity, all TE eigenvalues and all but the lowest TM eigenvalue are large and negative. Using 
(18), (19), (23), and (24), we have for the equivalent parallel-plane propagation constants 

kru^-j^ (31) 

^TE--jf- (32) 

These modes are all far below cutoff, and are attenuated in horizontal distances small com- 
pared with the vertical height of the ionosphere. The lowest TM mode, however, corresponds to 

ao=l (33) 
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^TMo — ^» (34) 

The introduction of an earth and ionosphere of reasonably high conductivity will not affect 
the eigenvalues drastically; ao will still be of the order of unity, although it will now be com- 
plex, and the remaining eigenvalues will be large and primarily real and negative. The low- 
est TM mode will still be the only one capable of propagation over significant distances. 
The following orthogonality relations are of interest in connection with the application 
of the radial wave functions to actual situations: 



X 



_„ € 



/: 



UaUa'dx=0 (a^T^a) (35) 



v^v^dx=0 (iSV/3). (36) 



The infinite limits represent integrations into fictitious ^^fiattened'' spaces above and below 
the surface of the earth. This artifice is permissible because of the high conductivity of both 
bounding media. 

An exact solution of either of the radial wave equations, or even an approximate solu- 
tion by any conventional method, is extremely difficult in any practical case. It is possible, 
however, to obtain a solution valid below the transition region and, in special cases, a solu- 
tion valid above that region. In the case of a perfectly conducting earth and finitely con- 
ducting ionosphere, the interior solutions are given by (27) and (28), as long as (ijOr/oo)<^l. 
If the effect of the earth is represented by a wave impedance Ze — this will be small compared 
to 770 — (27) and (28) become 

Ua= cos ^ll—ax—j ~ — (37) 

L 7?„V1— aJ 



Wl 

. z,Vi^" 



L ^0 J 



(38) 



The exterior solutions which can be obtained apply only where (o^r/oo) is large. One 
soluble case is that for which oor is given by a simple power law. If 

e=^l-jei{x-Xiy^-je,{x-Xi)'' (39) 

where n can have any real positive value, the solutions are 

Ua=A(x—Xiy^^-\^eH^i\ — j-;5 V^(^— ^i) (40) 

n+2 U^-T^ J 

v,= B{x-x,y''Hl_ [— ^^e(x-xA (41) 

n + 2 [_n-\-Z J 

The above expressions are restricted to the region where (a5;./co)^|Q:|. Another soluble case, 
which would appear to represent more nearly the true state of affairs, is the exponential 
variation of conductivity: 

e= 1 -iei6°^ - -jeie""'. (42) 

The corresponding exterior solutions are 

^« = ^V^fl^Vl+(4/a2)(a-l) (-Ve j (43) 

^^=5^(2)a)^lQV^) (44) 
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where the restriction is that (co^/co)^!. It should be noted that for the propagating mode 
(4/a^)(a— l)<d for the real ionosphere, and that therefore the order of the Hankel function 
becomes essentially unity for that mode. 

For many purposes, only the propagating mode is of interest. It happens that, in this 
exceptional case, a good method exists for joining the interior and exterior solutions, and cal- 
culating the eigenvalue. Because the permittivity varies very greatly over distances short 
compared with a wavelength, no WKB or quasi-WKB method is applicable in the transition 
region. However, the rapidity of this variation ensures that Ua^ will not change greatly Irom 
its value at x=0 until after the region of validity of the exterior solution is reached. For this 
reason, a good approximation to the solution in the interior and transition regions can be 
obtained by transforming (20) into an integral equation. If we let 

i^a(0)-l (45) 

then 

\e dx J^^.^ ' 7)0 
Note: Here and in most of what follows, a will be used to denote the lowest eigenvalue only. 
The first integration of (20) yields 

and a second integration leads to 

Ua{x) = l+j?^ {\{xyix'- {\{x') ['' Ua{x'')cWdx'+a {\{x') P' ^-^^^ rfx ' ' (Zo: ' . (48) 

^0 Jo Jo Jo Jo Jo ^(-^ ) 

If x is not too large, an iterative process starting with 

Ua{x'')^l (49) 

will converge to give as accurate a value of Uot{x) as nuiy be desired. The eigenvalue may then 
be calculated by setting the value of {lleUa)(dUa/dx) obtained by means of (48) equal to that 
obtained from (40) or (43) at some height Xq, chosen as small as consistent with the validity 
of the exterior solution. 

At sufficiently low frequencies, and of course depending upon the actual shape of the 
conductivity profile, the first-order approximation in which Ua{x) is simply unity can be used, 
and (47) gives 

1 dua Ze , f" dx' ,_^. 

e dx -^ rjo Jo €(xO ^ 

Under these circumstances, one obtains 

«= r:r:jz (51) 



r^odx 
Jo e 



where the last term in the numerator is to be calculated from the exterior solution. The 
term in Zg can usually be neglected at these frequencies. Moreover, the integral in the 
deiioiniriator is of the order of Xq, and so a more convenient form of (51) is 

I / 1 ^^« \ 
^ ^ \eUa dxjx^xp ^ . . 
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It should be remembered that (52) represents only a &st-order solution, but that the range of 
appHcabiUt}^ of the technique can be greatly extended by repeated iterations of (48). 

Equation (52) gives a means for calculating the lowest eigenvalue if the conductivity 
profile is known. Actually, however, the problem at hand is exactly the reverse — -that of 
calculating the profile if the eigenvalue is known as a function of frequency. The latter infor- 
mation may be obtained experimentally in a number of ways, two of which are described in 
the following sections. If a reasonable approximation for the upper ionosphere profile is 
known, (52) permits the calculation of the integral in the denominator as a function of fre- 
quency. This information is, in a general sense, a transform of the desired information on 
the variation of e with height. 

The desired information must be obtained by means of an inversion of this transform. 
Although the transform in question is not among those which have been well investigated in 
the past, several avenues of approach are open. It is conceivable, of course, that a rigorous 
analytical method for performing the inversion may be found. It is more probable, however, 
that recourse will have to be had to some approximate analytical technique. If all else fails, 
the transform can always be inverted by purely numerical methods. 

5. Field of a Vertical Dipole 

The atmospheric noise observed at ELF and below is due primarily to the vertical com- 
ponent of current in lightning strokes occurring at relatively large distances from the observer, 
and can therefore be described in terms of the field of an infinitesimal vertical dipole located 
just above the surface of the earth. This field is composed of TM waves of all orders with 
the azimuthal index m equal to zero. If the dipole is located at ^=0, x^O we have 

\p{x,d)=Y^ a.P.(-cos d)ua{x). (53) 

It is shown in the appendix that 

n"T) 

4 - uldx 

Jo € 

where p represents the current moment of the source. Except in the immediate vicinity of 
the source, only the lowest-order wave is significant, and therefore 

«'-''"=77-idTLT;r7^J«-w (55) 



j: 



dx 



where a denotes the lowest eigenvalue. When m = 0, as in this case, the only field components 
present are the following : 

Er=jv<,^4' (56) 

H.=i'4. (58) 

The eigenvalue can, in principle, be determined by comparative measurements of these 
three field components in individual sferics. Actually, insofar as the determination of iono- 
spheric properties is concerned, it is the quantity (a—1) which is significant. A comparison 
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of Er and H^ is not suitable for several reasons. First, and most important, is the fact that 
(a— 1) enters into the ratio only in second order. Second is the fact that the ratio depends 
to some extent on 6 — that is, on the distance of the observation point from the source. Third 
is the fact that the horizontal magnetic field measured by a loop antenna will depend upon the 
direction of incidence of the wave, while the vertical electric field will not. Consequently, 
two loop antennas will be required to determine the actual magnitude of H^, 

Comparative measurements of Ee and Er would yield (a — 1) in first order, although the 
experiment would have to be carried out at an altitude sufficiently great that Ee were indeed 
measurable. However, knowledge of the source distance and direction of incidence would 
still be required. 

The best procedure would be to measure and compare Ee and H^. The complex ratio of 
these two field components is entirely independent of the source distance, so long as this distance 
is great enough for the evanescent modes to be negligible. Moreover, the responses of two 
suitably oriented antennas will both depend in exactly the same way on the direction of inci- 
dence, so that the response ratio will be independent of the location of the source. An addi- 
tional important advantage is that the ratio is entirely unaffected by the presence of scattering 
centers such as discontinuities in the ionosphere (again provided the higher-order modes in 
the reflections can be neglected). The field-coin j)onent ratio is given by 



tIqII^ eUcc dx 

If the measurements are made in an airplane at dimensionless height Xa, (37) is applicable, 
the first-order approximation to tn^\^jl—a:x is valid, and e can be taken as unity, so that one 
has (with Zg — 0) 

-^=jXa{cc-\). (60) 

It should be emphasized that the values of a obtained in tliis way are those corresponding 
to the instantaneous conductivity profile above the point of observation. This method is 
thei'efore potentially capable of detecting both temporal and geographical variations in the 
ionosphere. In regions where the magnetic field of the earth is steeply dipping, the results 
obtained are simply those corresponding to the profile of (co,.)eff [Wait, 1960bl, rather than co^. 

6. Representation of the Field in Terms of Cavity Resonances 

Since (a— 1) is a quite slowly varying quantity, the field component ratio shown in (60) 
depends primarily on Xa or, for a given altitude in feet, on the first power of frequency. The 
same is approximately true of the ratio Ee/Er. It is this latter ratio which limits the measur- 
ability of Ee, due to the effect of errors in aircraft attitude. Assuming an altitude of 40,000 
feet and an accuracy of 0.1 degree in the determination of aircraft attitude, the probable error 
in Ee reaches 10 percent at about 50 cycles. Fortunately, the range below this frequency 
can be effectively covered through studies of the cavity resonances. A great deal can be 
inferred from ground-level measurements of the vertical electric field as a function of frequency 
in this range. However, such data can yield information on ionospheric properties only as 
averaged over the earth. More complete information could be obtained from simultaneous 
measurement of the other field components, but the following analysis is based on the vertical 
field alone, as these are the only data [Balser and Wagner, 1960a and b] presently available. 

The theoretical basis for this method is the following expansion : 

P.(-cos ^) _1 ^ (2r^+l)P,(cos d) . 

The vertical component of electric field at the surface of the earth, due to a single source, can 
be found from (45), (55), and (61). It is given by 
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^-<v)=-^^ w)=.(«)*+.) £ |!^;S:i-?j (62) 

where 

C(a>)= P^ (68) 

4x/*o — dx 

Jo e 

The total vertical field observed at a point is the result of some distribution of such sources 
over the surface of the earth. It is instructive to examine the case of a uniform distribution 
of a large number of noncoherent sources. The mean-square value of electric field is given by 

\E\ 2=^ ( '" f" ErE* sin eddd(t> (64) 

47r Jo Jo 

where Np is the total number of sources. Using (62) and the orthogonality properties of the 
Legendre polynomials, one obtains 

igp=iv.iaco)i-i: /f^+;^|'^^;+_^ii; - (65) 

For the rapidly varying ionosphere, C{cc) is roughly inversely proportional to co. 

7. Interpretation of Observed Atmospheric Noise Spectrum 

It was pointed out in the preceding section that airborne measurements would lose accuracy 
at lower frequencies and would cease to be valid below some minimum frequency. In order 
to determine that frequency, it was necessary to find the approximate relative magnitudes of 
the horizontal and vertical electric fields to be expected at flight altitudes. The noise spectrum 
observed by Balser and Wagner in the range from 5 to 34 cycles, which shows a series of pro- 
nounced peaks and valleys, has therefore been interpreted according to the scheme described 
below. The figure of 50 cycles as the lowest useful frequency for airborne measurements is a 
result of that anal^^sis. 

The development leading to (65) applies to a distribution of noncoherent sources at a 
single frequency. Actualh^, the noise energy is distributed throughout the frequency spectrum. 
The power per unit frequency interval received by an electrically short vertical monopole will 
be given by 

dP^Fjcc)^ (2n+l)(a'+b') 

do, co2 ^,[a-n{n+l)Y+b' ^^^^ 

where a and b are the real and imaginary parts of j'(?^+ 1), respectively, and i^(co) is a real function 
very nearly proportional to the frequency spectrum of the square of the current moment, 
(p^), of the sources. 

The problem of using (66) to find a as a function of frequency is not completely determi- 
nate, since there is an infinite number of ways to choose the frequency dependence of a, b, and 
F^ all of which will yield the observed spectrum. However, it is inconceivable that a and b 
could be anything but smooth, monotonic functions of w, and unlikely that i^(co), which includes 
geographical and time averaging of the sources, should be other than smooth and monotonic. 
The problem thus becomes that of choosing smooth forms for a and b such that the form of F 
then required to reproduce the spectrum is also smooth. 

The first step in the procedure is to plot cx)^(dPjdo)) as a function of frequency from the 
experimental data. The peaks in this function occur because of the resonances of successive 
terms in the series of (66), and will be located approximately at those frequencies where a is 
equal to n{n-i-l). This fact is used to plot several points in the (a,co) plane. A smooth 
curve is drawn through these points, and is represented as nearly as possible by a power law, 
or a combination of power laws. In this way, an initial approximation to a(a;) is found. Next, 
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it is assumed that both b and F are constant in the immediate neighborhood of each resonant 
peak. This permits an approximate determination of b, by means of (66), from the shape of 
each peak. The points tluis found provide a first approximation to b(co), which is also repre- 
sented as a power law. Equation (66) is then used to calculate the initial approximation to 
F{(x)) over the entire frequency range, which in general shows some irregularities. 

From this stage onward, the process is one of ^'cut and try^^ to find the smoothest possible 
F((x)). It does not appear to be possible, with these data, to obtain a perfectly smooth i^(co) 
using this model. There is a tendency for small apparent peaks of i^(co) to occur near resonances 
of even n, and small valleys near resonances of odd n. This should not be too surprising, 
however, as the time-averaged distribution of lightning strokes is not really uniform, but 
rather is more highly concentrated in a belt centered on the geographical equator. Such a 
distribution would be expected to couple more strongly to the even -numbered cavity modes. 
Furthermore, the fixed latitude of the observation point relative to this belt of excitation will 
tend to favor some modes over others. Still another possible source of discrepancies in the 
analysis is the fact that no account has been taken of the diurnal variation of e. However, as 
will be pointed out in the following section, the diurnal variation probably does not seriously 
aftect the vaUdity of the analysis. 

The expressions obtahied for a and b represeiit, to first order, an average of day and night 
conditions. The more laborious part of the analysis has not yet been carried to the point 
where the results can be presented with complete confidence. However, a reasonably good 
fit of the data was obtained with the following expressions: 

a=11.4(0' (/<20^) 



a=ll 



b 






(/>20~) 



(67) 



The expression for a is believed to be accurate within a very few percent, while that for b may 
be in error by as much as 10 percent, over the range 6 to 34 cycles. 

The corresponding source function F(o)) can be represented very well by 



7< = const.X / 1 + 8.38 exp [-3.15 (^)jy 



(68) 



This function is shown hi figure 2, and the theoretical and experimental noise spectra are 
displayed in figure 3. 



Figure 2. Source function F versus frequency. 
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Theoretical and experimental noise 
spectra. 



Figure 4. Phase velocity and attenuation versus 
frequency. 
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Because b^<^a^, it is possible to write for the equivalent plane-parallel propagation constant 



^TM — 



K-«^.) 



From this, one has 



'='Hi^) 



V 2.69 //\ 
c -yja \20/ 



db/1000 km 



(69) 

(70) 
(71) 



where A is the attenuation constant and v is the phase velocity. The values of these 
quantities obtained with (67) are shown as functions of frequency in figure 4. Although 
there is no reason to suppose that the empirical equations (67) should hold outside the range of 
the original data, it is nevertheless interesting to extrapolate them to higher frequencies. The 
extrapolated attenuation at 100-- is 1.375 db/1000 km, a value which is intermediate between 
existing experimental values [Jean, 1961] for da^^ and night at that frequency. The predicted 
phase velocity remains nearly constant over the whole range, and at 100'^ is equal to 0.797 
times the velocity of light. 
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8. Effect of the Diurnal Variation 



It is well kiiown that the properties of the ionospliere during (hvvliglit hours diiFer (jiiite 
markedly from those at uiglit. The quc^stion arises as to how the observed cavity resonances 
are ait'ected by this variation. In the absence of any diurnal variation, each cavity mode is 
associated with a pair of integer indices m and n^ and corresponds to a particular complex 
resonant frequency Smn (where ,y = ^+.7w). The modes are frequency degenerate, as -s-^^^^ 
depends only on n. The introduction of an angular dependence of e removes this degeneracy, 
splitting apart the resonant frequencies of the same n but different m. One would expect this 
to result in a certain amount of blurring of the observed resonant peaks with a consequent 
lowering of the apparent Q and raising of the apparent h. 

To a first approximation, the diurnal variation can be represented by a discontinuity in e 
across a plane which divides the spherical coordinate system in half at d = Trl2. In each hemi- 
sphere, the electromagnetic field can be represented in terms of solutions of the form described 
in (2) through (7). The boundary conditions which must be satisfied along this plane are that 
Hrj He, 7/^, Er, eEe, and E4, be continuous. When the expressions for the field components are 
combined with the partial chiferential equations satisfied by tfie potentials immediately to each 
si(h' of llie discontinuity, and accoimt is taken of the fact that in theneighboi'hood of 0=7r/2 the 
sine of 6 is both stationary and equal to unity, we are led to the following system of four fmic- 
tions, the continuity of whicli is sufficient to ensure that the boundary conchtions are satisfied: 



Vo 









dx . . 1 dV 



dd'^'^'''\ dxdct> 



a.^+^^ 






dx 



(72) 
(73) 
(74) 
(75) 



The apj)hcation of tlie above con(htions to the proldem of a vertical (hi)ole radiating at an 
arbitJ'ary point hi one of the hemispheres is extremely complex. The configuration is as sliowii 
in figure 5, where the source is located at A, a pohit not symnu^rically placed relative to the 
discontinuity. The total field consists of (1) a prinuuy TM wave set up by the source hi its 
own hemisphere and describable b}^ means of a potential iI/q, (2) reflected waves in the same 
hemisphere, and (3) transmitted waves in tlie other hemisphere. 

The primary TM potential ipo is independent of azimuth in the coordinate system x, 6' , </>', 
which has O^ as its polar axis. Provided A is not too close to the discontinuity, (55) can be 
used to compute this potential, which can then be expressed in terms of x, 6, </>, coordinates 



HEMISPHERE I 



DISCONTINUITY *- 



HEMISPHERE 2 




Figure 5. Configuration for study of diurnal 
effect. 
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centered on OB. The potential does depend upon azimuth in this coordinate system. The 
reflected and transmitted waves, in general, contain both TM and TE components. The 
potentials in the hemisphere containing the source are of the form 



i^r- 



kp 



r 



W.-oAi)^^ 



rP^Q(— cos ^0 cos ^— sin ^o sin ^ cos 0) 



sm vqit 



] 



^1,^0 (^) 



+Z)S«r COS m(/>Fr (cos d)ui, ,(x) (76) 



Xi-=SS*M sin 7n(t>P';^ (cos d)vi,^{x) 

n m 

while those in the other hemisphere are given by 

4^2=^2^^ COS m^Pf (-COS d)u2,K{x) 

K m 

^2=2Z;^x' sin 7n(t>P^(-cos e)i\x{x). 



(77) 

(78) 
(79) 



The reflected wave has +cos 6 and the transmitted wave —cos 0, as no singularity must 
appear at either ^=0 or 6=Tr. The indices 1 and 2 refer to the two hemispheres; vq is the 
lowest eigenvalue in Hemisphere 1; the v and /z are eigenvalues in Hemisphere 1, and k and X 
are eigenvalues in Hemisphere 2; (0, 6o, 0) are the coordinates of A in the system of OB. 

The coefficients in (76) through (79) can, in principle, be evaluated by application of the 
boundary conditions, which yields a separate set of equations for each value of 7n: 

jy)olZa':^Pf\^)uiAx)-m^h^I^{Qi)v[,,{x)+j^^ 

V iL K 

/»27r 

+ ^Zlc?xf'x(0)^2,x(^)=^m(^) COS ^0 P..(-sin (9o cos 0) cos m06?0 (80) 

X Jo 



^0- 



^0- 



im^S(i-Pr(0)<,(x)+S6-P;?'(0)z^i,,(x)-im^l]c-Pr(0)<,(^) 

+Z:<^r^'(0)«'2,x(x)=Z);(:c) ^i-^° P:„(-sin eocos0) siii</, sin m<^^ (81) 
X «i Jo 

l]M(M+l)&?P;(0)i'i.,(a;)-Z;X(X+l)fZrPr(0)»2.x(x) = (82) 

-i-SK''+l)arPr(0)Mi,.(:r)+i- S'<('<+l)crPr(0)M2.«(a;) 

=B^{x) ''"^'''"^^^ f" P.,(-sin 6>o cos 0) cos m^<t> (83) 
€i Jo 

where primes denote differentiation with respect to the argument of a function, and where 



Dm{x)=[l—^br,^A 



OJJJLoP 



/'oo 

47r| {ul,j€i)dx 



'ui,,q{x) 
_sin voTT 



]■ 



(84) 



Equations (80) through (83) have not yet been studied at great length. It can be seen, how- 
ever, that severe difficulties exist. A sufficient number of linear equations for the coefficients 
can be obtained by utilizing the orthogonality properties of the u or the v, but all these equations 
still have an infuiite number of terms. 



476 



The case in which no source is present and the two hemispheres are not very different has 
been briefly considered in an approximate way. The right-hand sides of all the equations are 
then zero, and oscillations can occur only at certain complex natural resonant frequencies. 
Each of the corresponding resonant modes has associated with it some particular value of m — 
that is, there is no mixing in m. For each value of m, modes are possible having an index n 
which may take on any integer value from m to infinity. These indices correspond to those 
of the P^(cos d) which one has in the ahsence of the discontinuity in e. There is, of course, 
mixing of eigenf unctions in v, /x, /c, and X. 

The highly tentative conclusions reached can best be explained with reference to the 
diagram of figure 6, which shows a portion of the complex frequency phme. Point E represents 
the resonant frequency corresponding to n=\l and a perfectly conducting ionosplierc located 
at some finite height. There is a four-fold degeneracy, since the modes for m=0, 1, 2, 3 all 
have the same resonant frequency. Point F represents the complex resonant frecpiency for an 
ionosphere having 

2eie2 



€l + €2 



(85) 



The four-fold degeneracy is still present. Points 0^ TI, </, and K show the way in which the 
introduction of the discontinuity separates the resonances of the same n and different m. If 
the vertical and horizontal shifts from E to F are regarded as first order, then, very roughly 
speaking, the horizontal separation of J and K is of third order, while the vertical separation 
is of fourth order. Resonances having (m+7i) odd tend to be shifted upward and to the right, 
while the reverse is true when (m-f-Ti) is even. The shifts seem too small to hnpair the validity 
of the analysis of section 7, although this should not be regarded as definitely established. 



PERFECTLY CONDUCTING IONOSPHERE- 
( degenerate) 



AVERAGE IONOSPHERE 
(degenerate) 




Figure 0. Nondegeneracy produced by diurnal 
variation of e. 



n = 3 



9. Appendix 

In order to evaluate the coeflicients a,, in (53), the fields on a small vertical cone surrounding 
the clipole and extending to infinity will be considered. It can be shown that 



and therefore that 



T • 7~» / /^x 2 sm VTT T 
Lim /^;,(— cos 6)= m 6j 

d->0 ^ 



Lim \l/{x,d) --=- In 6^a^ sin vw Ua{x) . 

0^0 TT V 



(86) 

(87) 
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In this region the magnetic field is 



H4>=- ^^=^S^^ sin vTTUa {x). {^^) 

V Ou wfu J, 

On the other hand, the magnetic field must be given by 

where j) represents the moment of the source, or the product of current and vertical current 
path: 

p=IL (90) 

Equating the two expressions for H^, one obtains 

y^.g. sin VTT Ua{x)==-^ 5(x). (91) 

The earth will be considered to be perfectly conducting, and Ua(0) will be taken equal to 
unity. Then (91) can be multiplied by Ua'{x)le and integrated. Making use of the fact that 
€(0) — 1, the orthogonality of the Ua then yields (54). 
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